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Galactic Cosmic Rays (GCRs) are mainly protons confined in the galactic 
magnetic field to form an isotropic flux inside the galaxy. Before reaching the 
Earth orbit they enter the Heliosphere and undergo diffusion, convection, mag- 
netic drift and adiabatic energy loss. The result is a reduction of particles flux 
at low energy (below 10 GeV), called solar modulation. We realized a quasi 
time-dependent 2D Stochastic Simulation of Solar Modulation that is able to 
reproduce CR spectra once known the Local Interstellar Spectrum (LIS). We 
were able to estimate the different behaviors associated to the polarity depen- 
dence of the Heliospheric modulation for particles as well as for antiparticles. 
We show a good agreement with the antiproton/proton ratio measured by 
AMS-01, Pamela, BESS, Heat and Caprice and we performed a prediction for 
the AMS-02 Experiment. 
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1. Introduction 

The effect of the heliospheric structure on GCRs propagation can be re- 
produced by a two dimensional (radius and helio-colatitude) Stochastic 
model solving numerically the Parkers's equation 1 . If we do not take into 
account the effects of the Earth magnetosphere 2 modulated fluxes depends 
not only on the level of solar activity but also on particles charge sign and 
solar magnetic field polarity 3 . The study of the modulation of p/p ratios 
is particularly important, because it includes explicitly the combination of 
charge sign and polarity dependence. The Local Interstellar Spectra (LIS) 
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used as input of the code, both for protons and antiprotons, are taken by 
the Galprop model 4 . 

2. Stochastic Model, Parameters and Data Sets 

2.1. Main parameters of the Model 

The present code simulates the interactions of a GCR entering the helio- 
sphere which extends up to a fixed distance of about 100 AU from the 
Sun. 

One of the main parameters of the model is the diffusion tensor. The 
parallel diffusion coefficient is: K^=ki(3Kp(P)(BQ/3B), where P is the 
particle rigidity (usually expressed in GV), k\ is the diffusion parameter 
discussed in the next section and Kp « P. The perpendicular diffusion 
coefficient has two components, radial K± r , and polar K±g. We used the 
relation: K± r = (K±)oKn, where {K±)q = 0.05. We also considered K±g = 
K± r in the equatorial region, while we enhanced its value in the polar 
regions of the heliosphere 5 : K±$ = 10.ffj_ r . 

We used the tilt angle a of the heliospheric current sheet (HCS) as a 
parameter for the level of the solar activity 6 : the higher the value of a the 
lower the expected GCR flux, for both solar field polarities. Values of the 
tilt angle are computed using two different models: the usual model uses 
a line-of-sight boundary condition at the photosphere and includes a sig- 
nificant polar field correction; an alternative model uses a radial boundary 
condition at the photosphere, and requires no polar field correction. As 
suggested by Ferreira and Potgieter 7 , the classical model is used for periods 
of increasing solar activity (for example 2007-2012, AMS-01 data, AMS- 
02 data), while the new model fits better for periods of decreasing solar 
activity (for example 2000-2007, BESS data). 

The three drift components do not depend on external parameters, ex- 
cept the solar polarity 8 (A>0 for positive periods and A<0 for negative pe- 
riods 6 ). The general drift expression is locally unlimited for a quasi- isotropic 
distribution 9 ' 10 , therefore we limit all the drift components below (n/A)v, 
which is the spatially averaged maximum value. 

2.2. Data Sets 

We selected GCR proton and antiproton data from 5 different experi- 
ments in order to compare and tune model results: AMS-01 11 , Caprice 12 , 
BESS 13-15 , HEAT and Pamela 16 The first two experiments took data in 
a period of positive solar polarity, BESS in both solar polarities, and the 
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last two in A<0 period. The corresponding periods of measurements are: 
June 1998 (AMS); May 1998 (Caprice); from July 1997 to December 2004 
(BESS); June 2000 (HEAT); and from 2007 to 2008 (Pamela). Solar wind 
values for these periods have been obtained from omniweb 17 by 27 daily 
averages, while tilt angle values from the Wilcox Solar Laboratory 18 . We 
can estimate the diffusion coefficient from a long term study 19,20 of neu- 
tron monitor measurements and the Force Field model 21 (FFM) approach 
to the treatment of Heliosphere transport of GCRs, using <~ 40 years of 
data ending in 2004. 
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Fig. 1. Diffusion parameter Kq - rising and declining phases for both negative and 
positive solar magnetic-field polarities - as a function of the SSN value; the continuous 
lines are obtained from a fit of Kq with respect to SSN values. 



In FFM 21 [e.g., see also Section 4.1.2.3 of Leroy and Rancoita (2009)], 
Gleeson and Axford (1968) assumed that i) modulation effects can be ex- 
pressed with a spherically symmetric modulated number density U of GCRs 
- the so-called differential density - with kinetic energy between T and 
T + dT, ii) the diffusion coefficient at the time t is given by a separable 
function of r (the radial distance from the Sun) and P (the particle rigidity 
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in GV): 

^r,t)=0fci(r,t)A>(P,t), 

with (3 = v/c , v the particle velocity, c the speed of light, Kp(P,t) w P 
for particle rigidities above w 1 GV and iii) the modulation occurs in a 
steady-state condition, i.e., the relaxation time of the distribution is short 
with respect to the solar cycle duration so that 

dU n 

They derived that the differential intensity at a radial distance r is given 
by the expression 



J(r,E t ,t) = J(r tm ,E t + * p ) 



P 2 - m r 2 c 4 



(1) 



.(^ + *p) 

where J(r tm ,Pt + 3>p) is the undisturbed intensity beyond the solar wind 
termination located at a radial distance r tm from the Sun; E t is the total en- 
ergy of the particle with rest mass m r and, finally, 3> p is the so-called force- 
field energy loss 21 ' 23 . When modulation is small 21 ' 23 - i.e., $ p -C m r c 2 ,T -, 
they determined that 

ZeP 
K P (P,t) 

where Ze is the particle charge and </> B (r,t) is the so-called modulation 
strength (or modulation parameter). Assuming that v w (the solar wind 
speed) and k\ are almost constant, <j> s (r,t) - usually expressed in units of 
GV (or MV) - reduces to 

, , v v (t) (r tm - r) 
0s(r '* j 3fei (t) ' (J 

from which one gets 

30s (r,i) 

i.e., fci is linearly dependent on (r tm — r). In the FFM, the diffusion coeffi- 
cient fC(r,t) is scalar quantity and, as a consequence, does not account for 
effects related to the charge sign of the transported particles. <j> a (r, t) is in- 
dependent of the species of GCR particles [e.g., see discussion at page 1014 
of Gleeson and Axford (1968) or Equation (1) of Usoskin and collaborators 
(2005)]. The values of the modulation strengths [0 S ('"Earth)] were monthly 
determined for the time period 20 from 1951 up to 2004 using measurements 
of neutron monitors (i.e., located at rEarth = 1 AU); while those for the 
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solar wind speeds are available on the web 24 . It has to be remarked that 
k\ depends on the value of the solar wind termination located at a radial 
distance r tm related, in turn, also to the solar wind speed [e.g., see Chapter 
7 of Meyer- Vernet (2007) and Section 4.1.2.2 of Leroy and Rancoita (2009)]. 
However, because the present simulation code assumes a fixed solar wind 
termination at 100 AU to calculate the modulated differential intensities at 
?"Earth, one has to derive from the diffusion parameter k\ that (K ) for an 
effective heliosphere with a radial extension of 100 AU (see Sect. 2.1). Thus, 
using Eq. (3) one can obtain 



99 AU 

K « fc 1 7 = 99 AU 

(rtm ~ fEarth) 



(4) 



_3r/> s (r E arth 

where 99 AU is the distance of the Earth from the border of the effective 
heliosphere as defined in the current simulation code. In Fig. 1, the diffusion 
parameter K Q - obtained from Eq. (4) - is shown as a function of the 
corresponding Smooted Sunspot Number (SSN) value. 




The Kq data had to be subdivided in four sets, i.e., rising and declining 
phases for both negative and positive solar magnetic-field polarities. For 
each set, the data i) could be fitted with a relationship - indicated in in 
Fi. 1 - between Kq and Smoothed Sunspot Numbers SSN 26 values and ii) 



December 23, 2010 7:40 



WSPC - Proceedings Trim Size: 9in x 6in grandi'll-nov-2010 



6 




■ I 




^10" 5 r 




■ PAMELA Adriani et. Al. 2010 - 



0.1 



1 



10 



100 



Energy (GeV) 



Fig. 3. Simulated p/p ratio as a function of the kinetic energy in GeV at 1 AU in 
comparison with experimental data: Pamela (2007-2008). 

exhibited a Gaussian distribution of difference of K values from the corre- 
sponding fitted values. The RMSs of the Gaussian distributions were found 
to be w 0.1339, 0.1254, 0.1040, 0.1213 for the phases rising with A < 0, de- 
clining with A < 0, rising with A > 0, declining with A > 0, respectively. 
In this way - i.e., from the relationships found -, we can use the estimated 
SSN values to obtain the diffusion parameter Kq at times beyond 2004. 
In practice, this procedure allows one to extend the ~ 40 years period by 
exploiting the linear relationship between the fitted Kq values and the SSN 
values (one of the main parameters related to the solar activity). Thus, 
we introduced in our code a Gaussian random variation of Kq with RMSs 
corresponding to those found for each subset of data. Results of the sim- 
ulation with and without the Gaussian variation are consistent inside the 
uncertainties of the code 1 . 

2.3. Dynamic Parameters 

Our code simulates the interactions of a GCR entering the heliosphere from 
its outer limit, the helio-pause, located - as already mentioned - approxi- 
mately at 100 AU 27 , and moving inwards to the Earth located at 1 AU. We 
evaluated the time t sw needed to SW to expand from the outer corona up 
to the helio-pause. Considering an average speed of 400 km/s it takes nearly 
14 months. While the time interval r ev of the stochastic evolution of a quasi 
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particle inside the heliosphere from 100 AU down to 1 AU is ~ 1 month 
at 200 MeV and few days at 10 GeV. This scenario, where r ev < t sw and 
t sw >> 1 month, indicates that we can not use fixed parameters (monthly 
averages) to describe the conditions of heliosphere in the modulation pro- 
cess. In fact at 100 AU, where particles are injected, the conditions of the 
solar activity are similar to the conditions present at the Earth roughly 14 
months before. Therefore we consider r ev negligible with respect to t sw and 
divide the heliosphere in 14 regions, as a function of the radius. For each 
region we evaluated 28 Kq, a and V sw , in relation to the expansion velocity, 
in a dynamic way. In the future the time spent by a GCR particle inside the 
heliosphere, as a function of the stochastic path and of the particle energy, 
will be also taken into account. 
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Fig. 4. Prediction of modulated p/p ratio as a function of the kinetic energy in GeV at 
1 AU for AMS-02 experiment (January 2011). 



2.4. Antiproton/ Proton: Comparison with Data and 
Prediction for AMS-02 

We performed the simulations using dynamic values of Ko, a and V sw . 
Results are shown in Fig. 2 and 3. Simulated fluxes with dynamic values 
show a very good agreement with measured data, within the quoted error 
bars. This happens both in periods with A>0, in comparison with BESS, 
and in periods with A<0, in comparison with Pamela. This means that our 
dynamic description of the Heliosphere improves the understanding of the 
complex processes occurring inside the Solar Cavity. 
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The periodic behavior of the heliosphere allows us to predict, with a 
certain level of precision, the parameters needed for a simulation related 
to a time in the near future. In order to get these data we considered 
the prediction of SSN from IPS (Ionospheric Prediction Service) of the 
Australian Bureau of Meteorology 29 . 

We concentrate our simulations on the AMS-02 mission 30 ' 31 that will 
be installed on the ISS in February 2011, and, in particular, at a time 
approaching the solar maximum: January 2012. We show in Fig. 4 the 
predictions of GCR modulation for the antiproton/proton ratio. 

3. Conclusions 

We built a 2D stochastic Monte Carlo code for particles propagation across 
the heliosphere. Present model takes into account drift effects and shows 
quantitatively a good agreement with measured values, both for positive 
and negative periods and for different particles and charge sign. This is 
relevant because particles with opposite charge sign undergo a different 
solar modulation 3 . We compared our simulations with antiproton/proton 
ratios measured by BESS and PAMELA. We used dynamic parameters 
values (K , a and V sw ) for the related periods, in order to reproduce the 
propagation of incoming GCR through magnetic disturbances carried by the 
outgoing solar wind. The dynamic description of the heliosphere and the 
forward approach seem to reproduce better the real physical propagation 
of GCR in the solar cavity. In order to have a more sophysticatcd model we 
need to introduce a dependence on the particle time spent in the heliosphere 
and a larger statistics of measured data during negative solar field periods, 
as AMS-02 will provide in the next years. Recent measurements 16 have 
pointed out the needs to reach a high level of accuracy in the modulation 
of the fluxes, in relation to the charge sign of the particles and the solar 
field polarity 32 . This aspect will be even more crucial in the next generation 
of experiments 30,31 . 
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